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I. 



INTRODUCTION 



Advances in the material sciences area have resulted in 
metals with high strength to weight ratios. The application 
of such material can result in structures with low levels of 
damping. Exposure of these structures to high levels of 
excitation creates a need for damping techniques which are 
both inexpensive and effective. Introduction of a 
distributed viscoelastic material with an incorporated 
constraining layer has gained widespread acceptance as a 
means of reducing vibration response. Many finite element 
codes have been developed which incorporate frequency 
prediction methods for use in design applications, reducing 
the need for experimental analysis. However, accurate 
behavioral predictions of applied viscoelastic material have 
remained costly, due to the frequency and temperature 
dependence in viscoelastic materials. 

The response of the simplest structures can be 
represented by continuous classical methods but the finite 
element approach remains the method of choice when complex 
physical systems must be analyzed. The use of a finite 
element method to find the system response requires 
modifications to capture the environmental dependencies of 
the applied viscoelastic layer. One standard method of 
conducting a frequency analysis of dampened structures is 
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through repeated application of the frequency response 
analysis over incremented steps in frequency. A second 
approach, using modal strain energy approximations in 
conjunction with the finite element method, has proved 
effective in predicting excitation responses of structures 
with a constrained viscoelastic layer. This paper 
investigates the two prediction methods described above as 
applied to plate structures. MSC/NASTRAN was chosen for the 
analyses because of its widespread availability and previous 
use in this area. The prediction methods were applied to 
two plate structures, one with the damping material having 
edge exposure to the environment and a similar plate with 
the damping material protected from the environment and 
possible corrosion, through the use of an edge weld. The 
edge-welded plate was evaluated before and after welding to 
determine possible damage to the viscoelastic layer. The 
prediction methods were evaluated for ease of application 
and prediction accuracy as compared to the experimental 
data. 
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II. THEORY AND BACKGROUND 



Continuous structures in the physical world possess 
distributed characteristics of mass, damping and stiffness. 
[Ref. l:pp. 140-141] Structures with such characteristics 
are said to contain an infinite number of degrees of 
freedom. The classical vibration analysis of such a 
continuous system involves the construction of a discrete 
mathematical model with a finite number of parameters to 
approximate the overall system response. Construction of a 
model of N discrete elements allows the discretized 
continuous system to be represented by the coupled matrix 
equations: 

[M] { q (t) } + [C]<4(t)} + [ K ] (q(t) } = {Q(t)> (2.1) 

Linear transformation and normalization of the undamped 
portion of Equation (2.1) and assumption of damping 
approximated by an appropriate diagonal matrix gives the 
uncoupled equations: 

(a(t)} + [ 2 ?Wn ]{a(t)} + [Wn ] < a ( t ) > = U(t)> (2.2) 

<q(t) >=[$]{ a(t) } (2.3) 
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where : 



{q(t) } = displacement vector 

[<j)] = mode shape matrix 

[ oo 2 ] = vector of natural frequencies 
n 

C = damping factor 

{l(t)} = transformed matrix of applied loads. 

Application of modal analysis through developed finite 
element codes have resulted in readily obtainable solutions 
to Equations (2.2) and (2.3) for complicated, lightly damped 
structures. [Ref. 2: Secs. 2. 3, 4. 3. 2; Ref. 3:p. 2] 

Recently, viscoelastic materials have been added to 

structures to improve their overall damping characteristics. 
However, the addition of these viscoelastic materials to 
design structures introduce both frequency and temperature 
dependent variables that must be incorporated into the 
finite element analysis. [Ref. 4:p. 2] 

Viscoelastic damping has been observed in many polymeric 
materials and has been directly attributed to the relaxation 
and recovery of long carbon atom chains of the polymer after 
deforming. [Ref. 5:p. 67] Because of their effect on 

molecular motion, temperature and frequency influence these 
"relaxation and recovery" properties. Thus, a complete 
understanding of the temperature and frequency dependent 
behavior of the material is essential . 
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Moderate changes in temperature or large changes in 
frequency result in definite predictable alterations in the 
viscoelastic material properties of shear modulus of 
elasticity and loss factor. Actual quantitative changes in 
the damping properties vary substantially with the chosen 
material. Qualitative changes are categorized into three 
descriptive regions; the glassy regions, transition region, 
and the rubbery region as shown in Figure 2.1. In a given 
viscoelastic material, all three regions are typically 
contained in two to three decades of frequency. [Ref. 5:p. 
91] 

Damping at relatively low temperatures for the given 
material is in the glassy region and exhibits a low loss 
factor and high shear modulus of elasticity. As the 
material temperature is increased the shear modulus begins 
to drop slowly as the loss factor rises. Sufficient 
temperature increases allow the damping material to enter 
the transition region, characterized by a rapidly decreasing 
shear modulus value. The loss factor reaches a peak value 
at the maximum rate of change in shear modulus value. 
Further increases in temperature result in the material 
entering the rubbery region with a corresponding decrease in 
loss factor and shear modulus values. While the above 
discussion is descriptive of varying properties with 
temperature increases, it also applies to frequency 
decreases. Examination of several materials shows a 25 
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degree fahrenheit increase being roughly equivalent to a one 
decade frequency decrease. For example, a 25° F increased 
temperature is roughly equivalent to a 10 kHz to 1 kHz or a 
1000 Hz-100 Hz drop in frequency. [Ref. 6:pp. A1-A21] 

The high sensitivity of viscoelastic materials causes 
the matrix values of Equation (2.2) to change with only mild 
variations in excitation and environmental conditions. This 
results in finite element solutions that are only applicable 
over a narrow range of frequency or temperature. A solution 
to Equations (2.2) and (2.3) which reflects varying matrix 
values may thus be obtained by the segmentation of the 
equations into incremental modal frequency response finite 
element analyses steps. This composite approach gives good 
results but in practical applications may become 
prohibitively expensive and time consuming, especially when 
several design alternatives are under consideration. 

An application of the modal strain energy method (MSE) 
by Johnson and Kienholz [Ref. 7;pp. 5-12] is suggested as an 
alternative solution to Equations (2.2) and (2.3) by 
imposing frequency dependency on the frequency response 
through post analysis data manipulation. In the MSE 
approach, the MS C/ NAS TRAN code is used to construct a finite 
element model with constant, midrange values of loss factor 
and shear modulus of elasticity. NASTRAN modal analysis is 
applied and the mode shapes and natural frequencies are 
extracted and inserted into Equations (2.2) and (2.3). The 
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fraction of elastic strain energy attributed to the 
viscoelastic core at each natural frequency is also obtained 
from the modal analysis and the damping material loss 
factors interpolated at each natural frequency. The matrix 

V V 

[2£ton] of Equation (2.2) is then replaced with [ri(^r-)Wn3 

V T 

when Equations (2.2) and (2.3) are modified using the 
interpolated loss factors, the frequency response solution 
becomes frequency dependent. The MSE method thus gives 
accurate solutions over a large frequency range while 
remaining cost effective. 
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Shear Modulus. G. or Loss Factor 




Temperature 



Figure 2.1 



Variation of Shear Modulus and Loss 
Factor with Temperature [Ref. 5:p. 90] 
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III. MATERIAL AND METHODS 



A. MODEL DESCRIPTION 

Three different plate structures were analyzed during 
this study in order to provide a basis for the comparison 
with the finite element results. These plate structures 

consisted of a simple plate, a milled plate with an 
unrestricted constraining layer and finally a milled plate 
with a restricted constraining layer. The simple plate was 
a 14 x 30 inch aluminum rectangular plate, 3/8 inch thick 
(Figure 3.1). Dimensions were chosen to give approximately 
15 resonant frequencies under 1000 Hz. The milled base 

plate and its dimensions are shown in Figure 3.2. All 

plates were machined from standard 6061 T6 aluminum. The 

use of a damped plate structure in a corrosive environment 
requires design modifications to protect the viscoelastic 
material . The welded plate structure of Figure 3.2 is such 
a design using the same material and overall thickness as 
the simple plate. This milled plate was subjected to 

vibration analysis before and after welding to detect 
reduced damping due to the welding process. 

B. DESIGN OF VISCOELASTIC AND CONSTRAINING LAYER 

THICKNESSES 

To approximate system loss factor, the plate structure 
was represented by a three layer beam system and subjected 
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to the Bernoulli-Euler type analysis described by M.L. Drake 
and developed by Ross, Kerwin, and Ungar [Ref. 8:pp. 1-6; 
Ref. 9]. The viscoelastic material and constraining layers 
were varied incrementally in the Bernoulli-Euler equations 
(3.1) to (3.3) below, using the FORTRAN program listed in 
Appendix A. 




2 



+ E_H_ (HL.-D) 



'3 3 ' 31 




- -?> + * E 3 H 3 H 31> 



(3.2) 





H 



H 1 +2H 2 +H 3 



■; and g = ~ 

E 3 H 3 H 2 K 



(3.3) 



31 



2 



where ; 



E = Young's modulus of elasticity. 



G = shear modulus. 



I = second moment of area. 



H = thickness. 



r / Beam Approximation Lengthy 2 A , Beam Approximation Widths 2 1 ^ 
K = Plate Length ' ' Plate Width 



and 
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Subscript l refers to the base structure. 

Subscript 2 refers to the damping layer. 

Subscript 3 refers to the outer layer. 

No subscript refers to the composite system. 

The maximum system loss factor and the temperature at which 
it occurred were found for each step. The results of 70 
incremental thickness changes were plotted in Figure 3.3. 

The University of Dayton Research Institute, in 
suggesting damping design procedures, recommends a system 
loss factor of 0.1 for most applications. [Ref. 10] 
Because Equations (3.1) -(3. 3) are only an approximation, the 
desired system loss factor was increased to 0.3 in this 
analysis. From this figure and based on a loss factor of 
. 3 , the values of the viscoelastic damping layer and 
constraining layer were chosen to be .0625 and .25 inches 
respectively . 

C. VISCOELASTIC MATERIAL 

The damping material applied to all three models tested 
was 3M SJ-201 ISD-112 viscoelastic compound with a thickness 
of 60 mil, supplied with suitable adhesive already applied 
to both sides by the manufacturer, 3M, St. Paul, Minnesota. 
During application of the viscoelastic and constraining 
layers, care was taken to allow air to escape from between 
the layers. Entrapped air can result in severe bond defects 



11 



giving lower natural frequencies and reduced system loss 
factor. [Ref. ll:pp. 8-9] 

ISD-112 was chosen for its excellent loss factor and 
shear modulus characteristics in the zero to 100° F 
temperature range for frequencies below 2,000 Hz. This 
trait is exhibited by the broad loss factor curve of Figure 

3.4 and is highly desirable in many applications. Figure 

3.5 shows a University of Dayton nomogram of the 
experimental determined shear modulus and loss factor 
properties for ISD-112. To determine shear modulus and loss 
factor values, locate the material temperature on the top of 
Figure 3 . 5 and the frequency of interest on the right 
margin. The diagonal line of temperature and horizontal 
line of frequency intersect on the vertical line of reduced 
frequency. Reduced frequency is the frequency value of 
interest multiplied by the temperature of interest and 
divided by the reference temperature of Figure 3.5. Reduced 
frequency and reference temperature are used to convert the 
input values of temperature and frequency to the conditions 
used to plot the loss factor and shear modulus curves. 
[Ref. 12:pp. 47-50] The intersection of reduced frequency 
and the plotted curves give the horizontal line of loss 
factor and shear modulus shown in the left margin. While 
nomograms are useful in comparing characteristics when 
selecting a damping material, precise values are difficult 
to read on the log-log scale and the actual equations. of the 
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nomogram curves were used to give accurate values of loss 
factor and shear modulus in the finite element codes. The 
equations for ISD-112 are given below and used in the 
program listed in Appendix A. [Ref. 9:pp. 6-14] 

Log(FR) = log(F) - 12* (T-TZERO) / (525*T-TZER0) (3.4) 

Log (G) = Log (ML) + (2*Log 10 (MROM/ML) ) / (1+FROM/FR) N (3.5) 

Log (ETA) = Log ( ETAFRO ) + A*B+D* ( 1-SQRT ( 1+A 2 ) ) / (C/2 ) (3.6) 

where the ISD-112 viscoelastic material nomogram constants 
are: 



TZERO 


= 131 


FROM 


= 400,000 


MROM 


= 2059.5 


N 


= .27 


ML 


= 28.427 


ETAFRO 


= 1.25 


SL 


= 0.27 


SH 


= -0.35 


FROL 


= 150,000 


C 

and: 


= 1.3 


FR 


= reduced frequency 


F 


= frequency 
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T = temperature 

G = shear modulus 

ETA = loss factor 

A = (Log (FR) - Log ( FROL) )/C 

B = SL+SH 

D = SL-SH 

D. NASTRAN FINITE ELEMENT MODEL 

Computer analysis of the sandwich structure was 
performed using the MSC/NASTRAN finite element program. As 
shown in Figures 3. 6-3. 8 the base plate and constraining 
layer were constructed of QUAD4 plate elements with 3/16 and 
1/8 inch offsets, respectively. The more complicated HEXA8 
solid elements used to simulate the sandwiched viscoelastic 
layer were required to accurately represent the strain 
energy due to shearing of the viscoelastic core and to 
account for the coupling between stretching and bending 
deformations. [Ref. 7:p. 8; Ref. 13:Sec. 2.22,4.3] The 

support bulk data card was used to give free-free end 
conditions to the plates. After constructing the finite 
element models, two different analyses were performed to 
predict frequency response amplitude. 

1 . Modal Strain Energy (MSE) Method 

Rearranging Equations (2.2) and (2.3) and solving 
for the response amplitude gives the following equation: 
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2 



x(ou) = 



? tl- (o 1 / Mn ) 2 }W E[A , 1 L )/ Mn 2 



V.. 



1 (1 - (w/wj 2 2 + ^n i (^) i (w/oo n ^ i )} 2 



n 

N 

+ I 



{ "i l §i l “\i lll *E,A,i 11 / "h 2 



2 \ 1/2 



1 {1 - (w/w ) 2 } 2 + (n. (-^) . (oj/oj .)}“ 

n X Vrp X -li/X 



V, 



where : 



X(oj) = amplitude of response as a function of 
frequency 

= ith natural radian frequency 

i = m °de shape value for the excitation grid 
point at the ith natural frequency 

< 1 >r i = mode shape value for the response grid 
point at the ith natural frequency 

nj_ = iterated viscoelastic material loss factor 

(V V /V T ) i = fraction of elastic strain in the visco- 
elastic material at each model 

L = magnitude of the load applied at the 
excitation grid point. 



The MSE method utilizes the MSC/NASTRAN method to find the 
constant values on the right side of Equation (3.7). The 
response amplitude is then calculated external to the finite 
element code. The mode shapes, natural frequencies and 
fractions of elastic strain energy in Equation (3.7) were 
determined using a single MSC/NASTRAN normal mode analysis 
for each plate. The damping material properties input to 
the MSC/NASTRAN normal analysis must be constant and were 
determined using the nomogram Equations (3. 4) -(3. 6) for a 
temperature of 60° F and the mid-range frequency of 500 Hz. 
The output of the MSC/NASTRAN normal analysis gave the 
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natural frequency, mode shapes, and the fractions of strain 
energy in the damping material at each associated natural 
frequency. The simple plate MSC/NASTRAN normal analysis 
program is listed in Appendix B. 

The loss factor of Equation (3.7) is sensitive to 
the natural frequency and several MSC/NASTRAN mode analyses 
must be applied to compensate for errors in the derived 
natural frequencies. Johnson and Kienholz describe an 
iterative method to approximate this loss factor frequency 
sensitivity. [Ref. 7:pp. 27-31] The iterative method was 

applied by first plotting the damping materials shear 
modulus as a function of frequency using the nomogram 
Equations (3.4)-(3.6). The damping material shear modulus 
of each plate was then varied through several MSC/NASTRAN 
modal analyses and the resulting change in each natural 
frequency was plotted on the same curve. The intersection 
of the curves gives an approximation of the natural 
frequency and shear modulus for each mode. Figure 3.9 
illustrates the iterate technique for two of the simple 
plate modes. The loss factor value was then found from this 
approximate shear modulus using the nomogram Equations 
(3.4)-(3.6). When this iterative method was first applied 
by Johnson and Kienholz, they noted the results differed 
from those observed in frequency response methods. Modes 
below the mid-range frequency show too much damping and 
those above show too little. [Ref. 7:p. 12]. Johnson and 
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Keinholz suggested multiplying the loss factor obtained in 
the described iterative method by (G^/G re f ) V 2 , where G is 
the viscoelastic material shear modulus at the ith natural 
frequency and G re f is the core shear modulus at the mid- 
range frequency. Once the matrices of natural frequency, 
mode shape, strain energy fraction and loss factor were 
calculated for each plate, Equation (3.7) was solved and the 
frequency response amplitude was plotted as a function of 
frequency. 

2 . Frequency Response Method 

Use of the NMS C / NAS TRAN model frequency response 
analysis to directly calculate response amplitude within the 
finite element program requires the viscoelastic material 
properties to be input as constants throughout the frequency 
range analyzed. This restriction limits the modal frequency 
response analysis to frequency spans over which the visco- 
elastic materials shear modulus and loss factor have little 
change. Inspection of the nomogram in Figure 3.5 indicated 
a frequency span of no greater than 100 Hz would give little 
change in the viscoelastic material properties. In order to 
give the best accuracy at the peak response values, the 
natural frequencies of each plate were estimated by 
conducting a MSC/NASTRAN modal analysis using the mid-range 
viscoelastic material properties. The viscoelastic material 
shear modulus and loss factor were then obtained for the 
plate structures at each of these natural frequencies. 
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Modal frequency response analysis was then performed over 
100 Hz spans centered at each natural frequency. The output 
magnitude versus frequency curves for each resonant 
frequency were compiled and plotted as composite curves. 
The modal frequency response analysis in MSC/NASTRAN for the 
simple plate is provided in Appendix B. 

E. VIBRATION GENERATOR TESTING ARRANGEMENT 

Experimentally determined vibration response analysis 
utilized a Wilcoxon research F3 vibration generator as the 
input force device. The F3 vibration generator was attached 
to the base plate on all applications. An HP 3562A Dynamic 
Signal Analyzer provided the swept sine driving signal which 
was amplified by a Wilcoxon research power amplifier (Model 
PA7C) . A force transducer was located before the attachment 
port at the base of the vibration generator to measure the 
input excitation force. The force signal was amplified by a 
charge amplifier and fed as the input signal to the HP 3562A 
Dynamic Signal Analyzer. An Endevco accelerometer was 
placed at a symmetrically opposite point on each plate and 
coupled with a signal conditioner to provide the output 
signal to the HP 3562A Dynamic Signal Analyzer. Figures 
3. 6-3. 8 show the locations of the vibration generator and 
accelerometer. The HP 3562A Dynamic Signal Analyzer 
converted the signals to digitized data and computed the 
frequency response, coherence and power spectra. 
Descriptive points were taken from the HP 3562A Dynamic 
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Signal Analyzer display of frequency response and plotted on 
the same graph as the finite element predictions. 
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Figure 3.1 Dimensions of Simple Plate with Applied 
Viscoelastic and Constraining Layers 
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Figure 3.2 Dimensions of Milled Plate with Applied 
Viscoelastic and Constraining Layers 
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Figure 3.4 Loss Factor Variation of ISD-112 with Varying Frequency and Temperature 
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Figure 3.5 University of Dayton' Nomogram for 3M ISD-112 Damping Material 





Figure 3.6 MSC/NASTRAN Finite Element 
Model for Simple Plate 
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Figure 3.7 MSC/NASTRAN Finite Element Model for 

Milled Plate with Free Constraining Layer 
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Figure 3.8 MSC/NASTRAN Finite Element Model for Milled 
Plate with Welded Constraining Layer 
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IV. RESULTS 



A series of response amplitude versus frequency plots 
were constructed for each plate. The first set of figures 
(Figures. 4. 1-4. 3) plots the experimental data for each 
plate before and after the application of the viscoelastic 
damping and constraining layers. After completing the two 
prediction methods, the results were plotted (Figures 4.4- 
4.6) with the experimental analysis for the damped plates. 

A. SIMPLE PLATE 

The experimental frequency response plot of the simple 
plate (Figure 4.1) shows a large increase in damping when 
the viscoelastic and constraining layers are applied. All 
natural frequencies shift lower in frequency and in some 
cases amplitude peaks are combined into a single peak. The 
system loss factor for the first mode increased from 2.2% to 
34% and the third mode increased from 1.6% to 14% after the 
application of the viscoelastic material and constraining 
layer. 

Comparison of the two simple plate prediction methods 
with each other and the experimental results show several 
differences (Figure 4.4). The MSE curve resonant 
frequencies differ from the NASTRAN composite curve, with 
the MSE resonant frequencies being slightly higher than the 
NASTRAN composite curve at frequencies below 500 Hz and 
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slightly lower at frequencies above 500 Hz. The resonant 
frequencies of the experimental data are lower than either 
prediction method throughout the frequency range inspected. 
The experimental data response peak amplitudes are higher 
than either of the prediction curves. The response peak 
amplitudes of the MSE method have approximately the same 
value as the NASTRAN composite curve with the poorest 
matching above 600 Hz and the best matching at the first 
peak. The experimental data show a loss factor of 34% for 
the first mode while both prediction methods indicate a 45% 
loss factor. The third mode has a 14% loss factor for the 
experimental data and a 20% loss factor for the prediction 
methods . 

B. MILLED PLATE, FREE CONSTRAINING LAYER 

The experimental frequency response plot shows a large 
decrease in peak response amplitudes after the viscoelastic 
damping treatment (Figure 4.2). In several cases, multiple 
amplitude peaks are again combined into a single peak. The 
system loss factor in the 300 Hz range was 1% prior to 
damping. Application of the damping material increases the 
system loss factor to 4% in this range. 

The prediction methods were plotted for the finite 
element model in Figure 4.5 with the two prediction methods 
showing good correlation. At frequencies below 500 Hz both 
prediction methods show similar peak response amplitudes but 
at differing frequencies with the NASTRAN composite method 
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having consistently lower peak frequencies. Above 500 Hz 
both prediction methods again show similar amplitudes but 
the MSE method showed lower peak frequencies. The 
prediction curves are poorly matched above 900 Hz. The 
experimental data show similar resonant frequencies but the 
frequency peaks are sharper than either prediction method 
over the observed frequency range. The experimental system 
loss factor at the 3 00 Hz peak is 4% while both prediction 
methods indicate a 6% loss factor. 

C. MILLED PLATE WITH WELDED CONSTRAINING LAYER 

After the experimental analysis of the unwelded milled 
plate, the plate constraining layer was welded to the base 
plate while maintaining the majority of the plate surfaces 
below 200° F during the welding process. After welding, the 
aluminum constraining layer and base plate were visibly 
warped and the plate resonated loudly when struck. The 
experimental analysis of the welded plate confirmed a 
significant decrease in damping (Figure 4.3) . The 
experimental loss factors of the welded plate at all natural 
frequencies were indistinguishable from the experimental 
loss factors for the milled plate prior to the application 
of the viscoelastic layer. AFter the weld was cut, the 
plate once again exhibited high damping (Figure 4.7). The 
plot of the two finite element predictions (Figure 4.6) 
yielded curves matched for frequencies below 700 Hz and 
curves beginning to disagree above 700 Hz. 
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Figure 4.1 Simple Plate, Experimental Data Before and After 
Application of Viscoelastic Damping Treatment 
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Figure 4.2 Milled Plate with Free Constraining Layer, Experimental Data 

Before and After Application of Viscoelastic Damping Treatment 




u 

w 



34 



Figure 4.3 Welded Plate with Welded Constraining Layer, Experimental Data 
Before and After Application of Viscoelastic Damping Treatment 
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Figure 4.4 Simple Damped Plate, Experimental Data and Two Predictions Method Curves 
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Figure 4.5 Milled Damped Plate with Free Constraining Layer, 

Experimental Data and Two Predictions Method Curves 
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Figure 4.6 Milled Damped Plate with Welded Constraining 
Laver. Two Predictions Method Curves 
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Figure 4.7 Milled Damped Plate Experimental Data Before 

Welding Constraining Layer and After Cutting Weld 



V. CONCLUSIONS 



Figure 4.1 shows the simple plate vibration response is 
reduced significantly after the 3M ISD-112 viscoelastic 
material and constraining layer are applied. Comparison of 
the experimental data and the two prediction methods in 
Figure 4.4, however, show that both the NASTRAN composite 
curve and the MSE method predict more damping than has 
occurred. Vibration energy absorbed by the viscoelastic 
layer is converted to heat and possibly increased the 
damping material temperature. No attempt to hold the plate 
temperature constant during the experimental analysis was 
made. Figure 3.9 shows that increasing the temperature of 
the viscoelastic material will decrease its loss factor. 
Another contributing factor to the lower actual damping may 
be the variation in the behavior of the viscoelastic 
material . 

In examining the MSE method (Equation (3.7)), it can be 
seen that the off-resonant frequency terms have an 
increasing contribution to the predicted response amplitude 
as the frequency increases. 

The natural frequencies of the MSE method are found 
using constant mid-range frequency viscoelastic material 
properties. This results in the MSE method amplitude peaks 
of Figure 4.4 shifting downward from the NASTRAN composite 
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curves in the high frequency range. An upward shift of the 
NASTRAN curves occurs in the low frequency range. The two 
prediction methods for the simple plate show more disagree- 
ment than for the other two plate models investigated. The 
MSE summation Equation (3.7) contained 18 modes for each 
plate analysis conducted. However the simple plate is the 
least stiff and thus 14 of these modes occur below 1000 Hz 
(78%) . The milled plate had 10 modes below 1000 Hz before 
welding (56%) and 8 after welding (44%) . The number of 
terms in the simple plate summation of Equation (3.8) should 
contain more terms to increase the accuracy of the MSE 
method prediction curve. 

Inspection of the experimental analysis in Figure 4.2 
again shows significant damping after application of the 
viscoelastic damping material and constraining layer. 
Figure 4.5 shows that both predicted damping levels were 
greater than the experimental data with the causes of the 
error being the same as for the simple plate. Examination 
of the two prediction curves shows the MSE method gives 
acceptable values of both amplitude and peak response 
frequency throughout the frequency range analyzed. The MSE 
prediction method of the milled plate with a free 
constraining layer again shows the MSE method response 
frequency shifts described in the simple plate MSE 
prediction. 
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After welding the constraining layer to the base plate, 
the structure was warped and exhibited no damping 
contributable to the viscoelastic material layer (Figure 
4.6). When the weld was cut, the plate warpage decreased 
significantly and the milled plate once again displayed a 
large loss factor (Figure 4.7). This indicated that the 
viscoelastic layer had separated from the plates as a result 
of the warpage caused by the edge welding process. Because 
of the loss of damping through deformation of the plates 
during welding, only the two prediction methods for the 
milled plate with a welded constraining layer are plotted in 
Figure 4.6. This plate was the stiff est plate analyzed, 
having eight natural frequencies in the frequency range 
studied vice 10 natural frequencies for the milled plate 
prior to welding and 14 for the simple plate. As a result, 
the response amplitude of the MSE prediction method is in 
agreement with the NASTRAN composite prediction method. The 
MSE prediction method again shows a slight upward shift of 
the low frequency range peaks and a downward shift of the 
high range peaks. 

Overall indications of the three plates show actual 
plate behavior is accurately estimated by both prediction 
methods. Errors occurring in the higher natural frequencies 
are due to an insufficient number of terms in the MSE 
summation Equation (3.4). Where inaccuracies occur, the MSE 
prediction method is found to be a good approximation of the 
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NASTRAN prediction method and should be an acceptable 
substitute for the NASTRAN composite prediction method. The 
simplicity of application coupled with the low computer time 
required make the MSE method the most attractive choice for 
engineering analysis which must consider both a wide range 
of environmental conditions and several design configurations. 
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VI. RECOMMENDATIONS 



The following are suggestions for further studies in 
this field: 

1. This study investigated prediction accuracies of the 

MSE and composite method at constant temperature. 
Accuracy over a temperature span should be 

investigated . 

2 . The results of the MSE method were heavily dependent 
on the MSC/NASTRAN normal mode analysis conducted 
with viscoelastic properties evaluated at 500 Hz. An 
investigation of the prediction sensitivity to shifts 
in evaluation frequency is recommended. 

3. The use of the HEXA8 solid finite element to model the 
viscoelastic layer gives results in an MSC/NASTRAN 
program which requires significant computer time to 
solve. Use of the viscoelastic beam element to model 
the viscoelastic layer have been suggested to reduce 
computer time requirements. [Ref. 14:pp. 352-354] 

4 . Failure of the welded plates viscoelastic material was 
due to the high temperature gradient causing warpage, 
yet failure to cool the surfaces of the plate would 
have resulted in destruction of the damping material. 
Investigation of welding techniques is required to 
give practicality to the milled plate design used in 
this paper. 

5. Application of several viscoelastic material and 
constraining layers have been used on structures in 
order to broaden the frequency range at which damping 
occurs. Applying the prediction methods investigated 
to such multiple damping layers should be researched. 

6. Viscoelastic material properties may vary from batch 
to batch of the same material. The construction of a 
nomogram using samples of the batch applied to the 
structure should greatly increase prediction 
accuracies . 
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APPENDIX A 



NOMOGRAM LOSS FACTOR AND SHEAR MODULUS CALCULATIONS 

As explained in Section III. A, it is impractical to use 
nomogram chart values of loss factor and shear modulus in 
analyzing damping characteristics of systems containing 
viscoelastic materials. All calculations of shear modulus 
and loss factor in this paper were accomplished through use 
of the FORTRAN program listed in Figures A.1-A.4. Although 
this program is not complicated, use of the equations 
without computer assistance is extremely tedious due to the 
large number of terms. Note that page 4 of the program 
(Figure A. 4) also gives the system loss factor based on the 
Bernoulli-Euler equations (3.2)— (3.3) • Equations (3.2) and 
(3.3) require Young's Modulus and the shear modulus values 
in order to be solved and thus it is most convenient to 
include these calculations here. 
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Figure A. 3 FORTRAN Program to Calculate the Viscoelastic Material Shear 
Modulus t Loss Factor and System Loss Factor (Page Three) 
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Figure A. 4 FORTRAN Program to Calculate the Viscoelastic Material Shear 
Modulus , Loss Factor and System Loss Factor (Page Four) 



APPENDIX B 



MSC/NASTRAN CODE FOR SIMPLE PLATE NORMAL MODE ANALYSIS 
AND MSC/NASTRAN CODE FOR SIMPLE PLATE FREQUENCY 

RESPONSE ANALYSIS 

While the MSC/NASTRAN code was used for all the finite 
element plates analyzed, listing each program would serve no 
useful purpose. For completeness, a normal mode analysis 
and a frequency response analysis are listed here for the 
simple plate structure (Figures B.1-B.10). The- MSC/NASTRAN 
programs for the milled plate with free constraining layer 
and milled plate with welded constraining layer vary only in 
the geometric organization of the grid points and finite 
elements . 
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ID DYNAMICS, FREQ MOD 

SOL 3 

TIME 240 

DIAG 3,9,13 

CEND 

TITLEsMODAL ANALYSIS 500 HZ PROPERTIES 

ECHO-UNSORT 

SPC - 1 0 

METHOD^ 1 

ESE=ALL 

DISP=ALL 

OUTPUT (PLOT) 

PLOTTER NASTRAN 
PAPER SIZE 14. BY 10. 

AXES=Z,X,Y 
SET 1 =ALL 

FIND, SCALE, SET 1, ORIGIN 1 

PLOT 

PLOT MODAU DEFORMATION 0, 1 , THRU, 15, SET 1, ORIGIN 1, SHAPE 
BEGIN BULK 
GRID, 1,0, 0,0,0 
-,*(7),=,=, *(3.75),= 

= 7 

GRID, 2, 0,2. 3333, 0,0 
=,*(7), =,=,*(3. 75),= 

= 7 

GRID, 3, 0,4. 6667,0,0 
— , * ( 7 ) ,— , — , *(3.75) ,— 

= 7 

GRID, 4, 0,7. ,0,0 
=,*(7), =,=,*(3.75),* 

= 7 

GRID, 5, 0,9. 3333, 0,0 
s ,*(7),=,=,*(3«75),s 
= 7 

GRID, 6, 0,11.6667,0,0 
*,*(7), =,=,*(3. 75),= 

= 7 

GRID, 7, 0,14. ,0,0 
=,*(7) ,=,=,*(3.75) ,= 

= 7 

GRID, 101, 0,0,0, .0625 
=,*(7), =,=,*(3.75),= 

= 7 



Figure B.l MSC/NASTRAN Normal Mode Analysis of 
the Simple Plate (Page One) 
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GRID, 102,0,2.3333,0, .0625 
=,*(7), =,=,*(3. 75),= 

= 7 

GRID, 103, 0,0. 5667,0, ,0625 
— ,*(7),=,=,*(3»75),= 

= 7 

GRID, 100, 0,7., 0, .0625 
=,*(7), =,=,*(3. 75),= 

= 7 

GRID, 105, 0,9. 3333,0,. 0625 
=,*(7), =,=,*(3. 75),= 

= 7 

GRID, 106,0, 1 1 .6667,0, .0625 
=,*(7), =,=,*(3. 75),= 

= 7 

GRID, 107, 0, ia.,0, .0625 



=,*(7), =,=,*(3.75),= 

= 7 

C3UAD0, 1 ,100, 1,2, 9, 8,, -.1075 
=,*(6),=,*(7),*(7),*(7),*(7),,= 

= 6 

C3UADa,2, 100, 2,3, 10,9, ,-. 1875 
=,*(6),=,*(7),*(7),*(7),*(7),,= 

= 6 

C3U ADO, 3, 10 0,3,9, 11, 10,, -.1875 
=,*(6),=,*(7),*(7),*(7),*(7),,= 

= 6 

C3UAD0,0, 100,4,5, 12, 1 1 , 1675 

=,*(6),s,*(7),*C7),*C7),*C7),,s 

= 6 

C3UAD9,5, 100, 5, 6, 13, 12,, -.1875 
=,*(6),s,»(7),*(7),*(7),*(7),,s 
= 6 

C3U ADO, 6, 100, 6, 7, 19, 13,, -.1875 
=»*(6),=,*(7),*(7),*(7),*(7),,= 

= 6 

C3U ADO, 1 01 , 101 , 101 , 102, 109, 10 8,, .125 
=,*(6),=,*(7),*(7),*(7),*C7),,= 

= 5 

C QU ADO, 102, 101 , 102, 103, 1 1 0 , 1 0 9 , , . 1 25 
=,*(6),=,*(7),*(7),*<7),*(7),,= 

= 6 



Figure B. 2 MSC/NASTRAN Normal Mode Analysis of 
the Simple Plate (Page Two) 
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C3UAD4, 103, 101 , 103, 104, 111, 110,, .125 
s ,*(6),=,*(7),*(7),*(7),*(7),,= 

-b 

C3UA04, 104, 101, 104, 105,1 12, 1 1 1, , .125 
= ,*(6),=,*(7),*(7),*(7),*(7),,s 
s 6 

C3UA04, 105, 10 1,1 05, 10 6, 11 3, 11 2,,. 125 
s,*(6),s,*(7),*(7),*(7),*(7),,s 
s 6 

CQUAD4, 106, 101, 106, 107, 114, 11 3,,. 125 
s,*(6),s,*(7),*(7),*(7),*(7),,s 
= 6 

CHEXA, 201 ,200, 1 ,2,9,8, 101, 102, +01 

*,*( 6 ),s,*( 7 ),*( 7 ),*( 7 ),*( 7 ),*( 7 ),*( 7 ),*(l) 

s ,*( 6 ),s,*( 7 ),*( 7 ),*( 7 ),*( 7 ),*( 7 ),*( 7 ),*(l) 

s,*< 6 ),s,*( 7 ),*( 7 ),*( 7 ),*( 7 ),*( 7 >,*( 7 ),*( 1 ) 

s,*( 6 ),s,»( 7 ),*( 7 ),*( 7 ),*( 7 ),*( 7 ),*( 7 ),*(l) 

s,*(6),=,*(7),*(7),*(7),*(7),*(7),*(7),*(l) 

=,*(6),=,*(7),*(7),*(7),*(7),*<7),*(7),*(1) 

s,*( 6 ),s,*( 7 ),*( 7 ),*( 7 ),*< 7 ),*( 7 ),*( 7 ),*(l) 

* 01 , 109,108 

* 1 , * 7 , *7 
*1 , *7, *7 
ft 1 , *7, *7 

* l , * 7 , * 7 
*1,*7,*7 
*1 ,*7,*7 

* 1 , *7 , * 7 

CHEX A, 202, 200, 2, 3, 10,9, 102, 103, f 09 

s,*( 6 ),s,*( 7 ),*( 7 ),*( 7 ),*( 7 ),*( 7 ),*( 7 ),*(l) 

s,*(b), 8 ,*( 7 ),*( 7 ),*( 7 ),*( 7 ),*( 7 ),*( 7 ),*(l) 

s,*( 6 ),s,*( 7 ),*( 7 ),*( 7 ),*C 7 ),*( 7 ),*( 7 ),*(l) 

=,*( 6 ),s,*( 7 ),*( 7 ),*( 7 ),*( 7 ),*( 7 ),*( 7 ),*(l) 

s,*( 6 ),s,*( 7 ),*( 7 ),*( 7 ),*( 7 ),*( 7 ),*( 7 ),*(l) 



s,*( 6 ),s,*( 7 ),*( 7 ),*( 7 ),*< 7 ),*( 7 ),*( 7 ),*( 1 ) 
s,*( 6 ),a,*( 7 ),*( 7 ),*( 7 ),*( 7 ),*( 7 ),*( 7 ),*(l) 
* 09 , 110,109 



ftl,*7,*7 

*1,*7,*7 

* 1 , *7 , * 7 

* 1 , *7 , *7 

* 1 , *7 , *7 



Figure B .3 MSC/NASTRAN Normal Mode Analysis of 
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* 1 , *7, *7 

CWEXA,203,200,3,4,11,10,103,104,+17 

s,*(6),s,*(7),*(7),*(7),*(7),*<7),*(7),*(1) 

=,*(6),=,*(7),*<7),»(7),*(7),*(7),*(7),*(1) 

= ,*(6),=,*(7),*(7),*(7),*C7),M7),*(7),*(1) 
=,*(6),=,*(7),*(7),*(7),*(7),*(7),*(7),*(1) 
=,*(6),=,*(7),*(7),*(7),*(7),*(7),*(7),*(1) 
=,*(6),=,*(7),*(7),*(7),*(7),*(7),*(7),*(1) 
=,*(b),:,*(7),*(7),*(7),*(7),*(7),*(7),*(l) 
♦17,111,110 

* 1 , *7 , *7 
*1,*7,*7 

* 1 , *7, *7 

* 1 , * 7 , * 7 
*1, *7, *7 

* 1 , *7, *7 

* 1 , * 7 , * 7 

CHEX A, 204, 200, 4, 5, 12, 1 1, 104, 105, +25 

=,a(6),=,*(7),*(7),*(7),*(7),*(7),*(7),*(1) 

s,*(6),s,*(7),*(7),*(7),*(7)'*(7),*(7),*(1) 

=,*(6),s,*(7),*(7),*(7),*(7),*(7),*(7)»*(l) 

=,*(6),=,*(7),*(7),*(7),*(7),*(7),*(7),*(1) 

*,*(6)»s,*(7),*(7),*(7),*(7),*C7),*(7),*(l) 

=,*(6),=,*(7),*<7),*(7),*(7),*(7),*(7),*(1) 

s,*(6),s,*(7),*(7),*(7),*(7),*(7),*(7),*(l) 

+25,112,111 

* 1 , * 7 , * 7 

* 1 , *7 , * 7 
* l , *7 , *7 
*1 , *7, *7 

* 1 , *7 , *7 
+ 1 , * 7 , * 7 

* 1 , *7 , *7 

CHEX A, 205, 200, 5, 6, 13, 12, 105, 106, +33 

=,*(6),=,*(7),*(7),*(7),*(7),*(7),*(7),*(1) 

s »*(6),s,*(7),*(7),*(7),*(7),*(7),*(7),*(l) 

=,*(6),s,*(7),*(7),*(7),*C7),*(7),*C7),*Ci) 

= ,*(6),=,*(7),*(7),*(7),*(7),*(7),*(7),*(1) 

=,*(6),=,*(7),*(7),*(7),*(7),*(7),*(7),*(1) 

=,*(6),=,*(7),*(7),*(7),*(7),*(7),*(7),*(1) 

=,*(6),=,*(7),*(7),*(7),*(7),*(7),*(7),*(l) 

+33,113,112 

* 1 , * 7 , *7 
*t,*7,*7 

* 1 , *7, *7 

* 1 , *7, *7 
*1 , *7, *7 

Figure B.4 MSC/NASTRAN Normal Mode Analysis of 
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*1 , * 7 , * 7 
* 1 , * 7 , * 7 

CHEXA,206,200,6,7, 14, 13, 106, 107, 
=,*(6),*,*(7),*(7),*(7),*(7),*(7),* (7 ),*(!) 



s,*(6),s,*(7),*(7),*(7),*C7),*(7),*(7),*Cl) 

s,*(6),s,*(7),*(7),*(7),*(7),*C7),*(7),*(l) 

s,*(b),s,*(7),*(7),*(7),*(7),*(7),*(7),*(l) 

=,*(6),=,*(7),*(7),*(7),*<7),a(7),*(7),*(1) 

=,*(6),=,*(7),*C7),*(7),*(7),*(7),*(7),*(1) 

=,*(6),=,*(7),*(7),*(7),*(7),*(7),*(7),*(1) 

♦41,114,113 

*1»*7,*7 

*1 ,*7, *7 

* 1 , *7 , *7 

* 1 , *7 , *7 

* 1 , *7, *7 
* i , *7, *7 
*1,*7,*7 

PSHELL, 100, 1 00, .375, 100 

P$HElL,101,100,.25,100 

PSOLID, 200,200,0 

MA71 , 100, 1 .*7, , .33, .0955 

MAT 1 ,200, ,2354.02, .49, .035,,, 1 .089 

PARAM, ASING, 1 

PARAM, AUTOSPC, YES 

SPC 1,10,6, 101, THRU, 163 

SPC1, 10,6, 1, THRU, 63 

SUPORT, 32, 12345 

EIGR, 1,MGIV, l.,1500.,,24,,,fZZ 
*ZZ,MASS 

PARAM, rtTMASS, .00259 
EMDDATA 



Figure B.5 MSC/NASTRAN Normal Mode Analysis of 
the Simple Plate (Page Five) 



54 



ID DYNAMICS, FREQ MOD 
SDL 30 
TIME 240 
CEND 

TITlEsFREQUENCY RESPONSE 500 HZ PROPERTIES 

ECHO=UNSORT 

SPC=10 

METHODS! 

SVECTORsALL 
DLOAD=2000 
FREQUENCY=2100 
SET 111=24 

DISPLACEMENT (PHASE )=1 1 1 
BEGIN BULK 
GRID, l, 0,0, 0,0 
=»*(7),=,=,*(3.75),= 

= 7 

GRID, 2, 0,2. 3333, 0,0 
s,*(7),=,=,*(3«75),= 

= 7 

GRID, 3, 0,4. 66b7, 0,0 
=,*(7), =,=,*(3. 75),= 

= 7 

GRID, 4, 0,7. ,0,0 
=,*(7), =,=,*(3. 75),= 

= 7 

GRID, 5, 0,9. 3333, 0,0 
=,*(7), =,=,*(3.75),= 

= 7 

GRID, 6,0, 1 1.6667,0,0 
=,*(7), =,=,*(3. 75),= 

= 7 

GRID, 7, 0,14. ,0,0 
=,*(7), =,=,*(3.75),= 

= 7 

GRID, 101,0,0,0, .0625 

*,*(7)., = , = ,*(3.75),* 

= 7 

GRID, 102,0,2.3333,0, .0625 
=,*(7), =,=,*(3.75),= 

= 7 

GRID, 103,0,4.6667,0, .0625 
=»*(7),=,=,*(3.75),= 

= 7 
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GRID, 104, 0, 7. ,0, .0625 
=,*(7), =,=,*(3.75),= 

= 7 

GRID# 105,0, 9.3333,0, ,0b25 
= , *(7),=,=, *(3,75),= 

= 7 

GRIO, 106,0, 11 .6667,0, .0625 
=,*(7 ),=,=, *(3.75),= 

= 7 

GRID, 107, 0,14.,0,. 0625 
=,*(7), =,=,*(3. 75),= 

= 7 

CQUAD4, 1 , 100, 1 ,2,9,8, ,-. 1875 
=,*(6),=,*(7),*(7),*(7),*(7),,= 

= 6 

C3UAD4,2, 100,2,3, 10, 9,, -.1875 

=,*(6),=,*(7),*(7),*(7),*(7),,= 

= 6 

C3UAD4,3,100,3,4,U,10,,-.l875 

=,*(6),=,*(7),*(7),*(7),*(7),,= 

= 6 

C3UA04, 4, 10 0,4,5, 12, 1 1,,-. 1875 
=,*(6),=,*(7),*(7),*(7), *(?),,= 

= 6 

C3UAD4,5, 10 0,5,6, 13, 12, 1 875 
=,*(6),=,*(7),*(7),*(7),*(7),,= 

= 6 

C3UAD4,6, 100, 6, 7, 14, 13,,-. 1875 
=,*(6),=,*(7),*(7),*(7),*(7),,= 

= 6 

C3IJAD4, 101, 101, 101, 102, 109,108, ,.125 
=,*(6),=,*(7),*(7),*(7),*(7),,= 

= 6 

CQUAD 4, 102, 101, 102, 103, 1 10, 109,, .125 
=,*(6),=,*(7),*(7),*(7),*(7),,= 

= 6 

C3UAD4, 103, 101 , 103, 104, 1 1 1, 1 10,, . 125 
=,*(6),=,*(7),*(7),*(7),*(7),,= 

= 6 

CQUAD4, 104, 101 , 104, 105, 1 12, 1 1 1, ,.125 
=,*(6),=,*(7),*(7),*(7),*(7),,= 

= 6 

C3UAD4, 105, 101, 105, 106, 1 1 3 , 1 1 2 , , . 1 25 
=,*(6),=,*(7),*(7),*(7),*(7),,= 

= 6 



Figure B.7 MSC/NASTRAN Frequency Response Analysis 
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CQUADfl, 106, 101,106, 107, 110,113, ,.125 
=,*(6 ),=,a(7),a(7),a(7),a(7),,s 
= 6 

CHEX A, 201, 200, 1,2,9,8,101, 102, tOt 

-,*(6),=,*(7),*(7),*(7),*(7),*(7),*(7),*(1) 

=,*(6),=,*(7),*(7),*(7),*(7),*(7),*(7),*(1) 

s,*(6),=,*(7),*(7),*(7),*(7),*(7),*(7),*(l) 

«,*(6),s,*(7),*(7),*(7),*(7),*(7),*C7),*C1) 

s,*(6),s,*(7),*(7),*(7),*(7),*(7),*(7),*(I) 

s ,*(6),s,*(7),*(7),*(7),*(7),*(7),*(7),*(l) 

s,*(6), = ,*(7),*(7),*(7),*(7),*(7),*(7),*U) 

♦01, 109, 108 

A 1 , *7, *7 

* 1 , * 7 f A 7 

* 1 , * 7 , * 7 
*1,a7,*7 

* 1 , *7 , a 7 

*1,*7,*7 

*1,*7,*7 

CHEXA,202,200,2, 3, 1 0,9, 1 02, 103, t09 

=,a(6),=,*(7),a(7),a(7),a(7),a(7),a(7),a(1) 

z,a(6),s,a(7),a(7),a(7),a(7),a(7),a(7),a(1) 

s ,a(6),s,a(7),a(7),a(7),a(7),*(7),a(7),a(1) 

s,A(6),s,A(7),A(7),*(7) f *(7),*(7),*(7) r *(i) 

=,*(6),s,a(7),a(7),*(7),a(7),a(7),a(7),a(1) 

s,*(6),=,a(7),a(7),a(7),a(7),a(7),a(7),a(1) 

s,a(6),z,a(7),a(7),a(7),a(7),a(7),a(7),a(1) 

♦09,110,109 

a 1 , a 7 , A7 

*1 , *7, *7 

* 1 , *7, *7 

*1 , *7, a7 

a 1 , * 7 , a 7 
A 1 , * 7 , A 7 

a1,*7,*7 

CHEXA,203,200,3,0,U,10,103,l04,tl7 

= ,A(6),s,A(7),A(7),*(7) f *(7) / *(7) ffc(7 )^ (l j 

=,a(6),=,a(7),a(7),a(7),a(7),a(7),a(7),a(1) 

=,*(6),=,*(7),AC7),A(7),A(7),A(7),A(7) f *(i> 

=,*(6),=,*(7),*(7),a(7),a(7),a(7),a(7),a(1) 

=,a(6),=,a(7),a(7),*(7),a(7),a(7),a(7),*(d 

=,*(6),=,A(7),A(7) f *(7) fA ( 7 ) # * (7)/ * (7)fjfc(1) 

♦17,111,110 
a 1 , *7 , * 7 
*1,*7,a7 
a 1 , a7 , a 7 
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*1,*7,*7 
*1 , *7, *7 

* 1 , * 7 , *7 

* 1 , *7, *7 

CH£XA,20a,200,«,5, 12, 1 1, 104, 105, +25 

=,*(6),=,*(7),*(7),*(7),*(7),*(7),*(7),*(1) 

s ,*(6),=,*(7),*(7),*(7),*(7),*(7),*(7),*(l) 

-,*(6),s,*(7),*(7),*(7),*(7),*(7),*(7),*(l) 

s,*(6),s,*(7),*(7),*(7),*(7),*(7)»*(7)r*(l) 

s,*(6),2,*(7),*(7),*(7),*t7)**(7)#*(7)»*(l) 

s,ft(6),s,*(7),*(7),*(7)»*(7),*(7),*(7),*(l) 

=,*(6),=,*(7),*(7),*(7)»*(7),*(7),*(7),*(1) 

♦25,112,111 

* 1 , *7, *7 
*1,*7,*7 

* 1 , * 7 , * 7 

* 1 , *7, *7 

* 1 , * 7, *7 

* 1 , *7 , *7 

* 1 , *7, *7 

C HE XA, 205, 200, 5, 6, 13, 12, 105, 106, ♦33 

=,*(6),=,*(7),*(7),*(7),*(7),*(7),*(7),*(1) 

s »*(6),s,*(7),*(7)#*(7),*(7)»*(7),*(7)#*(l) 

s,*(6),s,*(7),*(7),*(7),*(7),*(7),*(7),*(l) 

s,*(6),=,*(7),*(7),*(7),*(7),*(7),*(7),*(l) 

s#*(6),s,*(7),*(7),*(7),*(7),*(7),*(7),»(l) 

s,*(6),=,*(7),*(7),*(7),*(7),*(7),*(7),*(l) 

=,*(6),=,*<7),*(7),*(7),*(7>,*(7),*(7),*(1) 

♦35, 113,112 

*1**7, *7 

* 1 , *7, *7 
*1 , * 7 , * 7 
*1 , *7, *7 
*1,*7,*7 
*i,*7,*7 
*1 , *7, *7 

CHEX A, 206, 200, 6, 7,1 a, 1 3, 1 06, l 07 , ♦<* 1 

=,*(6),=,*C7),*(7),*(7),*(7),*(7),*C7),*(1) 

s,*(6),s,*(7),*C7),*(7),*(7),*(7),*(7),*(l) 

s,*(6),s,*(7),*(7),*(7),*(7),*(7),*(7),*(l) 

= ,*(6),=,*(7),*(7),*(7)»*(7),*(7),*(7),*(1) 

s ,*(6),s,*(7),*(7),*(7),*(7),*(7),*(7),*(l) 

=,*(6),=,*(7),*(7),*(7),*(7),*(7),*(7),*(1) 

s,*(6),=,*(7),*(7),*(7),*(7),*C7),*(7),*(1) 
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* 4i, ti4,i 13 

*lf*7/*7 

* 1 , *7 , *7 
*1,*7,*7 
*1,*7,*7 

* 1 , *7, *7 

* 1 , * 7 , *7 

* 1 , * 7 , *7 

DLOAD/2000, 1.0/ 1.0/2001 

RLOAD1, 2001, 2002, ,,2003 

DAREA, 2002, 40,3, 1 . 

TABLED1,2003,,,,,,,/+TA 

♦TA,0., 1 .,5000., 1 .0/ENDT 

FREQ1, 2100,450. ,5. ,20 

PSHELL, 100, 100, .375, 100 

PSHELL, 101 , 100, .25, 100 

PS0LI0, 200, 200,0 

MATl, 100, 1 .+7,, .33, .0955 

MAT1 ,200, ,2354.02, .49, .035,, , 1 .089 

PARAM, ASING, 1 

PARAM, AUTOSPC,YES 

SPCt, 10,6, 101, THRU, 163 

SPC1 , 10,6, 1 /THRU, 63 

SJPORT, 32, 12345 

EIGR , 1 , MGI V, 1 ., 1500. ,,20, , , tZZ 
tZZ,MASS 

PARAM, WTMASS, .00259 
ENDOATA 
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